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SUMMARY 

A study has been completed to determine the sensitivity of computed 

“"‘certainties in the Thermal Protection System 
vTPS) thermal model. Those parameters considered were; density, thermal 

insolation (RSI) 

uncei^aintr'^lOeT’'^”® thickness and emittance, and temperature measurement 
uncertainty. The assessment used a modified version of the computer program 

histori^ ^\hf orI^i!J“f Z calculate heating rates from temperature4ime 

• 1 original version of the program solves the "direct" one- 

d^ensional heating problem and this modified version of the program is set 
up to solve the ■'Inverse” problem. The modified program Is InteL™ Jor use 

modeinnralJererthI™“?“°H f" '“’"‘“I 

n!.! -f P ^ ^ thermal models were used in the study to determine the 

necessity for accurate knowledge of the thermal protection system's material 

-^eiri thermal properties the sensitivity (inaccuracies 

is very iL h convective heating rate by an altLed propertj) 

rpmii ^ ^ ^ great dependence on the front surface reradiative term 

requires accurate knowledge of emissivity and surface temperature fo^ar 
accurate heating rate determination. 


INTRODUCTION 

will Instrumentation (DFI) on the Space Shuttle Orbiter 

research opportunity to conduct flight aerothermodynamic 

esearch using the Orbiter as a flight research vehicle. Definition of the 
heating environment experienced by the shuttle's thermal protection systm 
throughout the entry trajectory is a significant part of this researX 
OrM^®^ number of thermocouples (approximately 225 ) distributed about the 
liltoVi "Z Z" measurement of surface tempLatuJe- ime 

t‘^^-J«‘:tory. These temperLure-tLr 

histories will be used in the calculation of convective heating rate hence 
defining the Orbiter’s entry heating environment. ® 

Determination of the convective heating rate to the surface of the TPS 
IvTllTu T transient, one-dimensional heat conduction 

!ni temperature history)! 

dlffpr is applied to a multi-layer, insulative TPS which is physically 

ifferent at each measurement point on the vehicle. Therefore, an existing 
lution method was sought which had multi-layer capability and was 
computationally efficient, and which provided the required solution accuracy. 


A search of various metho<)8 and available computer programs eliminated a 
number of possible methods (ref a 1 to 6) for various reasons. Some ware too 
expensive computationally (refs. 1 to 3), while others were not specific to 
the multi-layer insulatlve system problem (refs. 3 to 6) . The computer 
program described in NASA TM X-3370 (Pittman and Brinkley, ref. 7) was 
chosen for its low cost in computer time and its applicability to the shuttle s^^ 
layered TPS. Although this program provides a solution to the "direct problem, 
a simple modification to the finite difference equation at the heated surface 
allows for specification of the surface temperature-time history and solution 
for the convective heating rate. 

The accuracy of the calculated convective heating rate is dependent upon 
accurate knowledge of the location of the temperature measurement within the 
TPS, the thermal properties of the TPS materials, and the layer thicknesses. 
Using the Pittman and Brinkley program SINK in both the "direct" and Inverse 
formsi, the sensitivity of the calculated heating rate to uncertainties in 
thermal properties, layer thicknesses, and surface temperature has been 
determined. 


SYMBOLS 


specific heat, J/kg-K 


heat capacity of front surface heat sink (pC^t) , J/m -K 

total enthalpy at edge of boundary layer, J/kg 

local enthalpy of fluid at front surface temperature, J/kg 


‘*CALC 


thermal conductivity, W/m-K 

q as calculated by the inverse method, W/m 
c 


VrIVEN 


input to direct method, W/m 

2 

cold-wall convective heating rate, W/m 

2 

radiant heating rate to the surface, W/m 

4 2 

reradiated heating rate 1 ^1 


temperature, K 


temperature of first station of layer, K 


2 


t heat sink thlcknessas» m 

X coordinate normal to surface 

ot absorptance of front surface 

e . emlttance of front surface 

8,1 

0 Stefan-Boltzmann constant, 

T time, sec 


ANALYSIS 


Description of Method 


The original version (ref, 7) of the SINK program provides a one-d:^mensional 
analysis of the transient thermal response of multi-layer insulative systems. 

The surface heating rate, initial temperature distribution through the material, 
and material properties and layer thicknesses must be specified for solution. 

The program can handle multiple layers of material and air gaps between layers 
and at the back surface* The program can have both radiative and convective 
heating rate inputs. Its most applicable feature to the inverse problem Is 
that the temperature at any one station in the TPS material can be specified 
as a function time. With the front surface temperature specified (station 1), 
a simple modification to the front surface boundary condition allows 
calculation of convective heating rate. In finite difference form, the 
convective heating rates at the beginning and the end of each time step are 
assumed equal, as are the enthalpies. This assumption gives tlie res *lt shown 
below in equation (1). (Terms are not shown differenced here.) 
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For this calculation, q- is assumed zero. This is the only modification to the 
SINK program. 


Thermal Models 

Three different thermal models were used in the conduct of this study. The 
models provide mathematical representations of each of the three TPS types 
present on the Orbiter. The TPS types and thermal models are described below, 
and In figures 1 to 3. 
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High-Temperature Reusable Surface 

Insulation (HRSI) - Figure 1 

Consists of a Silica 

Model: 


tile (with a hlgh-emlttance. 

Coating 

4 Nodes 

black, water-proof coating) 

Tile 

86 Nodes 

bonded to a felt strain 

Adhesive 

2 Nodes 

isolator pad (SIP), which 

SIP 

2 Nodes 

is in turn bonded to the 

Adhesive 

2 Nodes 

vehicle structure. 

Structure 

2 Nodes 


Low-Temperature Reusable Surface Insulation (LRSI) - Figure 2 


Consists of a Silica tile 

Model : 


(with a high-reflectance, white, 

Coating 

A Nodes 

water-proof coating) bonded to 

Tile 

86 Nodes 

a felt SIP, which is in turn 

Adhesive 

2 Nodes 

bonded to the vehicle structure. 

SIP 

2 Nodes 


Adhesive 

2 Nodes 


Structure 

A Nodes 


Flexible Reusable Surface Insulation (FRSI) - Figure 3 


Consists of a blanket of 

Model : 


felt (with a high-reflectance. 

Coat ing 

A Nodes 

white, water-proof coating) bonded 

Felt 

86 Nodes 

to the vehicle structure. 

Adhesive 

6 Nodes 


Structure 

A Nodes 


Layer thickness dimensions shown In figures 1 to 3 are applicable to the 
specific body points modeled for the purpose of this studv. Body point 
locations are shown in figures 4 and 5. Heating rate and pressure Information 
was for trajectory 14 4 14, 1C. The Appendix gives thermal property information 
used in this study. 


Accuracy of Method 

In order to assess the accuracy of the inverse solution technique, the 
direct solution was applied to known heat ing-^rate histories for each of the 
selected body points to generate surface temperature-time histories, using nominal 
thermal models. These temperature-time histories va^re then used as Input to 
the Inverse solution, again with nominal thermal models. In order to generate 
computed heating rate histories. Comparison of the computed heating rate 
histories with the original given lieatlng-rate hlstorit^s provides a measure of 
the accuracy of the solution technique. This comparison is shown in figures 6 
to 9 for two of the body points. Solution error is not significant for 
trajectory times between 100 and 1500 seconds for any of the body points studied. 
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Sensitivity Analysis 

After the accuracy of the modified method was established, the method was 
used to parametrically determine the sensitivity of the calculated heating rate 
to uncertainties in thermal model parameters. The surface temperature-time 
histories generated using the direct solution with nominal thermal models were 
used as inputs to inverse solutions in which thermal model parameters were 
systematically varied. 

The sensitivity of the solution to an uncertainty in a thermal model 
parameter appears as an error in the computed convective heating rate when 
compared to the input heating rate of the direct solution. Table I gives a 
list of the parameter variations considered in the study. 


RESULTS AND DISCUSSION 

Although all variations listed in Table I were investigated, most had an 
insignificant effect upon the computed result; therefore, only those found to 
be significant will be discussed. The calculated convective heating rate for 
the HRSI system is affected by uncertainties in the surface emittance and the 
surface temperature. A 5-percent uncertainty in the surface emittance creates a 
5— percent error in the calculated convective heating rate. This is shown in 
figure 10. A 5-percent change in surface temperature (thermocouple error), 
figures 11 and 12, creates a 20-percent error in the convective heating rate 
calculation. The source of these errors is illustrated in figure 13 which shows 
a plot typical of the HRSI system with the convective heating rate and the 
reradiative term plotted versus time. The dominance of the reradlative term in 
the equation for q (aq. (1)) is evident from this plot. Consequently, a 
5-percent change in^emlttance causes a 5-percent change in the reradiative term, 
and a 5-percent change in surface temperature is magnified by the fourth power 
in the reradiative term. 

The LRSI and FRSI systems showed similar results (figs. 14 to 19) but 
additionally, significant errors in the calculated heating rate were noted from 
0 to 400 seconds and from 1200 to 1600 seconds in the trajectory, when the 
coating thicknesses were douoled (figs. 20 and 21). However, in the high 
heating phase of the trajectory, errors in the coating thickness will create 
errors in the calculation of less than 5 percent. No other parameter 
variations significantly affect the convective heating-rate calculation. 

A direct solution was obtained for the various RSI tile systems with 
temperature data from 0.076 cm (0.03 in) below the surface stored for use in the 
inverse calculations. These interior temperatures were used in the heating rate 
calculations as assumed surface temperatures to assess the errors resulting from 
a thermocouple which is not in thermal contact with the surface. The resulting 
errors in heating rate calculations are shown in figures 22 to 27. Results 
indicate that for similar types of tiles and similar depth of thermocouples, the 
solution gives similar error In heating rate for different heating environments. 
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This indicates that a correction night be applied to shuttle data based on 
type of tile and thermocouple depth for the various heating mvlronmenCs » if 
thermocouple depth is known* 

An investigation to determine the effect of uncertainty in the initial 
temperature distribution through the material on the calculated heating rate 
was also completed. Both high and low backface temperatures were Investigated 
with the temperature varying linearly through the material to the given surface 
temperature. A 55 K variation from front surface to back was imposed. The 
effect on the calculation for convective heating rate was insignificant. 


CONCLUDING REMARKS 

A study has been completed to determine the sensitivity of computed 
convective heating rates to uncertainties in the TPS thermal model. The 
assessment used a modified version of the SINK program described in NASA 
TM X-3370 to calculate heating rates from temperature-time histories. Both 
nominal thermal properties and alte-ed thermal properties were used. The 
modified program is Intended for use in thermocouple data reduction for shuttle 
flight data. The necessity for accurate knowledge of the thermal properties 
and layer thicknesses of the shuttle's thermal protection system was determined 
by the sensitivity study. Results of the study show very low sensitivity for 
many thermal properties; however, a great dependence on the front surface 
reradlative term was shown for accurate calculation of the convective heating 
rate. This dependence requires accurate knowledge of the surface emissivity 
and surface temperature for an accurate convective heating rate calculation. 

A slight dependence on the surface coating thickness was also determined for 
the LRSI and FRSI systems. Therefore, thermal properties of the shuttle's 
insulative system can vary from nominal properties (those specified by the 
manufacturer) somewhat without affecting the calculation of convective heating 
rate. However, due to the affect of the reradlative term in the calculation, 
accurate knowledge of the surface temperature and emissivity is desired for 
an accurate convective heating-rate calculation. 
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HRSI Tiles 


Appendix 

Material Properties 


C oating 

Thickness 

Density 

Emlttance 

Absorptance 


3.81 X lO"^ m 

1666 kg/m^ 

0.85 

0.85 


Temperature 

(°K) 


Thermal Conductivity 
(W/m-K) 


Specific Heat 
(J/kg-K) 


255.6 

533.3 

811.1 

1088.9 

1227.8 

1366.7 

1450.0 
1644.4 

1811.1 
1922.2 


0.843 

1.045 

1.218 

1.378 

1.449 

1.528 

1.565 

1.687 

1.870 

2.042 


794.96 

1004.16 

1192.44 

1317.96 

1380.72 

1443.48 

1476.95 

1569.00 

1631.76 

1631.76 


11. H-900 RSI 

Thickness 
Density 


7.94 X 10 ^ m 
144.17 kg/m3 


Temperature 

(°K) 

0 

Thermal Conductivity (W/m-K) 
Pressure (N/m2) 

10.05 100.5 lOlJ.O 10131.0 

101314.6 

Specific Heat 
(J/kg-K) 

255.6 

0.0130 

0.0130 

0.0173 

0.0317 

0.0433 

0.0476 

627.6 

394.4 

0.0160 

0.0160 

0.0216 

0.0390 

0.0547 

0.0600 

878.64 

533.3 

0.0216 

0.0216 

0.0290 

0.0478 

0.0700 

0.0750 

1054.37 

672.2 

0.0303 

0.0303 

0.0374 

0.0562 

0.0850 

0.0924 

1150.60 

811.0 

0.0403 

0.0403 

0.0476 

0.0680 

0.1040 

0.1140 

1205.00 

950.0 

0.0533 

0.0533 

0.0606 

0.0850 

0.1255 

0.1353 

1238.46 

1088.9 

0.0720 

0.0720 

0.0800 

0.1068 

0.1514 

0.1630 

1255.20 

1227.8 

0.0980 

0.0980 

0.1056 

0.1330 

0.1835 

0.1960 

1263.57 

1366.7 

0.1270 

0.1270 

0.1353 

0.1630 

0.2200 

0.2350 

1267.75 

1533. 3 

0.1670 

0.1670 

0.1765 

0.2008 

0.2700 

0.2900 

1267.75 


111. RTV-560 (Room Temperature Vulcanizing) 

Thickness 1.78 x 10~^ m 

Density 1409.7 Kg/m^ 

Thermal Conductivity 0.3115 W/m-K 

Specific Heat 1464.40 J/kg-K 
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IV. SIP 

Thickness 4.064 x 10 
Density 86.5 kg/m^ 


Temperature 

(°K) 


Thermal Conductivity (W/m-K) 
Pressure (N/m^) 



0 

10.05 

255.6 

0.0092 

0.0092 

311.1 

0.0098 

0.0098 

366.7 

0.0109 

0.0109 

422.2 

0.0126 

0.0126 

477.8 

0.0157 

0.0157 


100.5 1013.0 10131.0 
0.0190 0.0308 0.0343 
0.0215 0.0360 0.0407 
0.0234 0.0415 0.0472 
0.0263 0.0471 0.0550 
0.0291 0.0524 0.0642 


Specific Heat 

(J/kg-K) 

101314.6 

0.0355 794.96 

0.0422 1079.47 

0.0493 1439.30 

0.0571 1882.80 

0.0661 2405.80 


V. RTV-560 (same as III) 

VI. Aluminum 2219-T8XX 

Thickness 6.35 x 10“3 m 
Density 2803.4 kg/m^ 


Temperature 

(OK) 


Thermal Conductivity 
(W/m-K) 


Specific Heat 
(K/kg-K) 


255.6 

311.1 

366.7 

422.2 

477.8 

533.3 


119.42 

128.07 

135.00 

141.92 

146.60 

150.57 


899.56 

928.85 

953.95 

979.06 


LRSI Tiles 
I • Coating 


II. 


Thickness 

Density 

Emittance 

Absorptance 

Thermal Conductivity 

Specific Heat 

LI-900 RSI 


3.05 X 10'^ m 

1666.0 kg/m3 
0.80 
0.32 

same as HRSI Coating 


Thickness 0.014 -i 

All Thermal Prop. os same as HRSI 

III. RTV- 560 - same as HRSi RTV-560 


8 





IV. SI£ - same as HRSI SIP 

V. RTV-560 - same as HRSI RTV-560 

VI. Aluminum 2219-T8XX - same as HRSI item VI 

FRST Tiles 


I. Coating 


Thickness 

Density 

Emittance 

Ausorptance 

Thermal Conductivity 

Specific Heat 


1.76 X lO"^ m 
1553.86 kg/m3 
0.80 
0.32 

0.3115 W/m-K 
1A64.4 J/kg-K 


II. Felt 

Thickness 4.064 x 10 ^ m, 8.128 x 10 ^ m 
Density 86.5 kg/nj3 


Temperature Thermal Conductivity W/m-K 


(°K) 



Pressure (N/m^) 




0 

1.005 

10.05 

100.5 

1013.0 

10131.0 

101314.6 

255.6 

0.0138 

0.0138 

0.0182 

0.0242 

0.0300 

0.0343 

0.0356 

311.1 

0.0149 

0.0149 

0.0208 

0.0287 

0.035 

0.0412 

0.0433 

366.7 

0.0164 

0.0164 

0.0240 

0.0336 

0.0415 

0.0476 

0.0502 

422.2 

0.0176 

0.0176 

0.0270 

0.0384 

0.0476 

0.0560 

0.0580 

477.8 

0.0190 

0.0190 

0.0294 

0.0433 

0.0547 

0.0640 

0.0665 

588.9 

0.0225 

0.0225 

0.0360 

0.0545 

0.0704 

0.0822 

0.0846 

700.0 

0.0260 

0.0260 

0.0433 

0.0658 

0.0865 

0.1052 

0.1073 

311.1 

0.0303 

0.0303 

0.0520 

0.0800 

0.1064 

0.1340 

0.1376 

III. 

RTV-560 - same as HRSI RTV-560 




IV. 

Aluminum 221S 

1 

H 

00 

1 

same as 

HRSI item VI 




Specific Heat 
(J/kg-K) 


1305.41 
1388.90 

1401.64 

1443.42 
1506.24 
1589.92 
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Coating 3.81 x 10 " m 


Tile 7.9 X lO’^ m 



Figure 1.- HRS I thermal model (not drawn to scale). 


3.05 X 10 m 



Figure 2.- LRSI thermal model (not drawn to scale). 















1 

Felt 8.41 X 10'^ or 4.05 x lO’^tn 

- t 

Structure — 

6.35 X lO'^m 


Coating 
2.54 X lO’^m 

Adhesive 
1.78 X 10'^ m 



Figure 3.- FRSI thermal model (not drawn to scale). 
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Body point 1150 (HRSI) 
Body point 1650 (HRSI) 



Figure 5.- Body point locations used. Upper lower view. 




Figure 6. - Body poi 





Figure 7.- Body poi 
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Figure 14.- Body point 3629, LRSI, 5 percent change in emittance, error 
between direct and inverse methods. 
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Figure 20.- Body point 3629, 
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Figure 26.- 



Figure 27.- 



